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Abstract 

The maximal superintegrability of the intrinsic harmonic oscillator potential on N- 
dimensional spaces with constant curvature is revisited from the point of view of 
sZ(2)-Poisson coalgebra symmetry. It is shown how this algebraic approach leads to 
a straightforward definition of a new large family of quasi-maximally superintegrable 
perturbations of the intrinsic oscillator on such spaces. Moreover, the generalization 
of this construction to those TV-dimensional spaces with non-constant curvature that 
are endowed with sZ(2)-coalgebra symmetry is presented. As the first examples of the 
latter class of systems, both the oscillator potential on an iV-dimensional Darboux 
space as well as several families of its quasi-maximally superintegrable anharmonic 
perturbations are explicitly constructed. 



1 Introduction 



The Poisson-coalgebraic "dynamical" symmetry underlying all the superintegrable Hamil- 
tonian systems that we shall present in the sequel can be summarized as the following 
quite general result [1, 2]: Let (q, p) = ((gi, . . . , qw), (pi, . . . ,Pn)) be N pairs of canonical 
variables with respect to the Poisson bracket 

= £ (^-a£flfc). (L1) 



i=l 



x Based on the contribution presented at the "17th Conference on Nonlinear Evolution Equations and 
Dynamical Systems" NEEDS 2007. L'Ametlla de Mar, Spain, June 17-24, 2007. 
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and let us consider the three functions given by 

N N N 

i=l i=l i=l 

Then, given any smooth function 7i, the Hamiltonian 

i^)=W(q 2 ,p 2 ,q.p) (1.3) 

defines an iV-dimensional (iVD) classical superintegrable Hamiltonian with (2N — 3) 
functionally independent integrals of the motion that are explicitly given by 

m N 

c (m) = (q iPj -q jPi ) 2 , C\ m) = Yl (QiPj-QjPi) 2 , 

l<i<j N-m+l<i<j 

where m = 2,...,N and = C^y Furthermore, the sets of functions given by 

{H (N \C^} and {H^ N \ C (m) } (m = 2, . . . , N) define two sets of N integrals in involution. 
Proofs, technical details and further generalizations can be found in [1, 2] but, at this point, 
some remarks concerning the symmetry and superintegrability properties of are in 

order. 



• Remark 1. If the three functions (1.2) are written as J_, J+ and J3, respectively, 
by computing the Poisson bracket (1.1) among them we recover the Lie-Poisson 
commutation rules of sl(2,R): 

{J 3 ,J+} = 2J+, {J 3 ,J_} = -2J_, {J_,J+} = 4J 3 . (1.5) 

In other words, the functions (1.2) define a particular ND symplectic realization 
of st(2,R). Hence any Hamiltonian can be thought of as a smooth function 

defined on s((2, R) 

HW = W(J_, J+, J 3 ) = H(q 2 , p 2 , q • p), (1.6) 

and all the results here presented can be interpreted in the framework of sl(2,R)- 
Poisson dynamics. Properly speaking, is denned on a three-dimensional st(2, R)- 
subalgebra of the sl(2, R) (g) ■ ■ ■ (g> N ^ sl(2, R) Poisson algebra. We stress that integrable 
systems on the Euclidean space and endowed with the A^-particle st(2, R)-symmetry 
given by the representation (1.2), were already studied in [3]. 

• Remark 2. The "universal" integrals of motion (1.4) are derived from the Casimir 
function of the aforementioned s((2,R) Poisson algebra (see [1]), and are given as 
sums of the square of certain angular momentum components. In particular, since 
the functions Ljj = qiPj — qjPi with i < j and i,j = 1, . . . , N span an $o(N) Lie- 
Poisson algebra with Poisson brackets 

{Lij,Lik} = Ljk, {Lij,Ljk} = —Lit, {Lik,Ljk} = Lij, i < j < k, (1.7) 

the integrals (1.4) can be rewritten as 

TO N 

= E L h C {m) = £ L 2 ,. (1.8) 

l<i<j N-m+l<i<j 
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This, in turn, means that the sl(2, R) <8> • • • (8)^ sl(2, R) symmetry gives us the right 
prescription to get the appropriate subset of quadratic functions of the generators of 
so(iV) that Poisson-commute with the Hamiltonian (1.3) and are in involution. In 
this respect, note that the algebra sl(2,R) ® ■ ■ ■ (%> N ^ sl(2,R) has only 3iV generators 
and many of them do Poisson-commute, whilst so(iV) has N(N — l)/2 generators 
with many non-vanishing Poisson brackets among them. 

• Remark 3. It is well-known that the maximum number of functionally independent 
(and different from the Hamiltonian itself) integrals of the motion for an ND Hamil- 
tonian is (2iV — 2). In the case that all these integrals do exist, the sytem is called 
maximally superintegrable (MS). Since 

(1.3) has, by construction, (2N - 3) 
functionally independent integrals, we shall say that this is a quasi-maximally super- 
integrable (QMS) Hamiltonian. Nevertheless, for some specific choices of the function 
TC it will be possible to find the remaining integral (which is not provided by the 
above symmetry). In that case H will define a MS system. 

• Remark 4. The canonical variables (q, p) have a priori neither a given geometrical 
(physical) meaning, nor restricted (real/complex) values. 

Therefore, we can conclude that HW (1.6) comprises a large family of QMS Hamilto- 
nians; each particular system arises for a specific choice of the function H together with 
an "appropriate" geometrical interpretation of the canonical variables. 

1.1 Oscillators on the iVD Euclidean space 

In order to illustrate these ideas and also as the starting point for further developments, 
let us consider the iVD isotropic harmonic oscillator with angular frequency to. Such a 
system can easily be identified within the family (1.6) by simply setting 

n = \ J+ + \^ 2 J- = \v 2 + ^V, (i.9) 

for which q are Cartesian coordinates in the ND Euclidean space E^. This Hamiltonian 
is not only QMS, but the standard prototype of MS systems. In fact, the "remaining" 
constant of the motion can be taken as any of the N integrals 

li=p 2 i +oo 2 ql i=l,...,N, (1.10) 

since each Xj is functionally independent with respect to both the set (1.4) and H. The 
results above summarized allows for a straightforward superintegrable even-order anhar- 
monic oscillator perturbation given by [1]: 

11 00 11 °° 

k=l k=l 

which is QMS for any choice of the 5k parameters, since the Hamiltonian (1.11) is indeed a 
function of (1.2). We remark that once a single anharmonic contribution with parameter 
(5^ is added to the first harmonic term the maximal superintegrability is lost, but the 
resulting system (for any number of arbitrary <5fc's) always keeps the (27V — 3) integrals of 
motion (1.4). In particular, the latter are just the integrals of the motion for the radial 
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Garnier system [4, 5], which is recovered by taking to and 5\ as the only non- vanishing 
parameters. We also recall that the integrability properties of some quartic oscillators can 
be generalized to the Calogero-Moser systems defined with such nonlinear oscillators as 
external potentals (see [6] and references therein). 

In the following sections we present the superintegrable Hamiltonians defining an- 
harmonic oscillators on the ND sphere and hyperbolic spaces [7], as well as on an iVD 
Riemannian space of variable curvature. The latter space is an ND generalization of the 
2D Darboux surface of type III [8, 9], one of the four 2D spaces with non-constant cur- 
vature whose geodesic flows are MS. We stress that in the constant curvature cases the 
zero-curvature (flat) limit leads to the Euclidean nonlinear oscillator (1.11). 



2 Anharmonic oscillators on spaces of constant curvature 

In this section we consider the ND classical Riemannian spaces with constant sectional 
curvature re: the sphere 8^ (re > 0) and the hyperbolic space M N (re < 0). We recall 
that both of them can be embedded in a linear space IR^" 1 " 1 with ambient or Weierstrass 
coordinates (xq, x) = (xq, x±, . . . ,xn) subjected to the "sphere" constraint 

E : xl + kx 2 = 1. (2.1) 

The metric on the proper ND spaces reads [10]: 



ds 2 = - (dxo + redx 2 ) 



(2.2) 



where dx 2 = J^iLi dxf- 

In order to be able to apply the results described in section 1 to the geodesic flows 
and oscillator potentials on and M N (with the construction of the constant curvature 
counterpart of the QMS Hamiltonian on (1.11) in mind), we shall proceed as follows: 



• We interpret in a proper way the "abstract" canonical coordinates and momenta 
(q, p) as intrinsic quantities on each space; this step can be achieved through different 
projections from the ambient space M Ar+1 . The resulting metric in terms of q leads 
to the kinetic energy term of the Hamiltonian. 

• We deduce the form of the corresponding intrinsic harmonic oscillator potential on 
the spaces with constant curvature as a function of q; its flat limit (or contraction) 
re — ► has to give (1.9). Moreover, in our framework the well-known maximal 
super integrability of this potential on these spaces has to be explicitly proven by 
finding the additional integral of the motion through direct computation. 

• Finally, QMS anharmonic oscillator potentials can be obtained as a symmetry- 
preserving perturbation of the intrinsic oscillator on these spaces: in particular, 
we can consider the sum of all the powers of the above intrinsic oscillator potential 
in such a way that their flat limit re — > reduces to (1.11). In this way, the constant 
curvature analogues of the radial Garnier system will be obtained by considering the 
anharmonicity given by the square of the intrinsic oscillator. 
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In the following we apply the above steps by considering two types of phase spaces 
(q, p) in and M N coming from different projections from M JV+1 . Obviously, the two 
Hamiltonians so obtained are canonically equivalent through a change of coordinates de- 
spite its apparent disequivalence as objects defined on the sl(2, R) Poisson algebra. 



2.1 Stereographic projection: Poincare coordinates 

Let us consider the stereographic projection [7] from the ambient coordinates (xo,x) £ 
£ C R N+1 to the Poincare coordinates q € 1^ with pole (— 1, 0) € R N+1 : 

(-l,0)+A(l,q)€S. (2.3) 

Hence we obtain that 

X = TTkq*' - = A - 1 = ^?- x = A " = W < 24 > 

Therefore, the metric (2.2) in Poincare coordinates reads 

da 2 

dg2 = 4 n 2VT ( 2 - 5 ) 
(1 + Kq 2 ) 2 

And the associated geodesic flow has (up to a positive constant factor) the free Lagrangian 

t = 2otWf (2 - 6) 

The canonical momenta p are obtained through the usual Legendre transformation and 
read: 

Thus the ND kinetic energy on and M N is given in Hamiltonian form as a particular 
case of (1.6); namely 

T = I (1 + K J_) 2 J + = 1 (1 + K q 2 ) 2 p 2 . (2.8) 

As a consequence, this geodesic flow is (at least) QMS, with (1.4) being the explicit (and 
invariant) form of the integrals of the motion. However, we stress that in the context of 
Poincare coordinates the geometric interpretation of such integrals in terms of angular 
momentum components is lacking. 

The next point is to deduce the curved harmonic oscillator potential in terms of these 
Poincare coordinates q. We recall that the radial (geodesic polar) distance r from an 
arbitrary point to the origin in (k > 0) and M N (k < 0) along the geodesic joining 
both points is written (in ambient coordinates) as [10, 11] 

1 x 2 

-tan 2 (V«r) = — . (2.9) 

Consequently, the well-known intrinsic oscillator on constant curvature spaces (the so- 
called Higgs oscillator [12, 13]) is written in Poincare coordinates as 

2 
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In this way the full Higgs oscillator Hamiltonian reads 

« = \ (1 + K J.f J + + ^ IT -i rp = i (1 + «,»)' p> + \J ^±^. (2.1!) 

This system is known to be MS [12]. Therefore, the remaining functionally independent 
constant of the motion does exist and, therefore, it has to be found by direct methods. 
Such an additional integral can be shown to be any of the following N functions [1]: 

2 2 

X % = ( K (l-Kq 2 )+2K(q-p) g ,) 2 + - q \ i = l,...,N. (2.12) 

(L — KC[ ) 

Now, a natural perturbation of this Hamiltonian including anharmonic terms that 
preserve the QMS properties of the system would be 

« = \ a + u + ^ +|> ( <r^rp) M 

i i 2 00 2(k-\-l) 

= 1 (1 W) V + ^ + E 4 (1 .V,,^,, ■ ("3) 

Notice that when any <5fc 7^ this curved anharmonic oscillator is QMS (it always commutes 
with the integrals (1.4) due to its s[(2, R)-coalgebra symmetry) but not MS (at least with 
integrals depending quadratically on the momenta). In fact, the N = 2 restriction of this 
perturbed system does not appear in the classifications of MS systems on S 2 and H 2 given 
in [14, 15, 16]. Note also that the first perturbative term given by 5\ 7^ can be considered 
as the constant curvature generalization of the (radial) Gamier system. 



2.2 Central projection: Beltrami coordinates 

Now we consider the central projection from the ambient coordinates (xo,x) G S C R^" 1 " 1 
to the Beltrami ones q G WL N with pole (0, 0) G R N+1 : 

(0, 0) + fi (1, q) G S, (2.14) 

so that we find 

1 q 

/i=- 7 ==, x = ^, x = ^q = (2.15) 
V 1 + Kq 2 V 1 + ^q 2 



Then the metric (2.2) turns out to be 



(1 + /tq 2 )dq 2 - /c(q • dq) 2 
(1 + K q 2 ) 2 



The corresponding free Lagrangian is given by 

(1 + Kq 2 )q 2 - re(q • q^ 2 
2(1 + «q 2 ) 

which leads to the definition of the conjugate momenta p as: 

_ (1 + Kq 2 )q - re(q • q)q 
P " (1 + Kq 2 ) 2 • 



T = ^— — Z' H „„ 2 T HJ , (2-17) 



(2.18) 
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Hence the kinetic energy Hamiltonian describing geodesic motion reads 

T = \ (1 + kJ-) («/+ + «J|) = 1(1 + ^q 2 ) (p 2 + «(q • P) 2 ) • (2.19) 

By taking into account (2.9) we find that the expression of the curved oscillator on E> N 
and M N adopts in these Beltrami coordinates the following "Euclidean" form: 

U = ^V, (2.20) 

which yields the following expression for the complete curved oscillator Hamiltonian (again 
as a particular case of (1.6)): 

H = \ (1 + «J_) (J+ + kJI) + X -u? J_ = 1(1 + K q 2 ) (p 2 + «(q • p) 2 ) + l^ 2 q 2 . (2.21) 

The remaining constant of the motion for this MS Hamiltonian can be taken from any of 
the N functions [1] 

Z l = ( Pl + K{q-p)q^ 2 +cu 2 ql i=l,...,N. (2.22) 

And the explicit QMS anharmonic generalization of (2.21) is proposed to be 

1 1 °° 

H = - (1 + «J_) ( J+ + kJ 2 ) + -u; 2 J_ + <5 fe J^ +1 

k=i 

i 1 00 

= 1(1 + -q 2 ) ( P 2 + «(q • P ) 2 ) + l^ 2 q 2 + ^ . (2.23) 

fc=l 

Note that in this coordinates the curved Gamier term is given just by the q 4 perturbation. 



3 Oscillators on an 7VD space of non-constant curvature 

In arbitrary manifolds with non-constant curvature, kinetic- energy Hamiltonians can ex- 
hibit extremely complicated dynamics and are, in general, no longer integrable [17]. From 
a physical viewpoint, the caracterization of (super)integrable geodesic flows on curved 
(pseudo)-Riemannian manifolds in arbitrary dimensions is relevant for supergravity and su- 
perstring theories, and constitutes an active research field (see [18] and references therein). 
Let us consider the (spherically symmetric and conformally flat) iVD Riemannian manifold 
whose metric and geodesic Hamiltonian flow are given by 

d S 2 = (a + q 2 )dq 2 , T = (3.1) 

where the parameter a > 0. This is a space with non-constant curvature. Moreover, its 
scalar curvature is negative and given by 

R=-(N- ^ N ;^" N . (3.2) 

In the N = 2 case, this space is just one of the so-called 2D Darboux spaces: the 2- 
manifolds with non-constant curvature admitting two quadratic first integrals, so that its 
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geodesic motion is quadratically MS. There are only four types of such spaces [8], which 
are (from the integrability viewpoint) the closest ones to constant curvature ones, since 
their AD generalizations are the only spaces other than K N , M N and S N whose geodesic 
motion could be expected to be (quadratically) MS for any dimension. In fact, an AD 
spherically symmetric generalization of the four 2D Darboux spaces has been recently 
introduced in [19] and shown to be QMS. 

The integrability properties of the space (3.1) have been thoroughly studied in [20], 
where it has been shown that the Hamiltonian 

J+ + J 2 J_ p 2 9 q 2 / \ 

n= + - T — = ^-a+w -i-a, 3.3 

a + J_ a + q 2 a + q 2 

is a MS system with (2 A — 2) functionally independent quadratic first integrals given 
again by (1.4) and one of the following functions 

li= P 2 -{H-u 2 )q 2 l , i=l,...,N. (3.4) 

Note that this integral cannot be written as a function of the s((2, R) symmetry; moreover, 
the set {Xj : 1 < i < A} is also in involution. To the best of our knowledge, this 
Hamiltonian provides the first example of a Hamiltonian system on a Riemannian space 
of non-constant curvature which is MS in any dimension. 

Furthermore, a geometric analysis shows [20] that the potential 



J 2 q 2 

T~ = w 9 

a + </_ a + q z 



U = oj 2 — = to 2 (3.5) 



can be interpreted as the intrinsic harmonic oscillator on this curved space, that turns 
out to be MS, despite of the introduction of a non-constant curvature. We remark that, 
as expected, for A = 2 this model is listed in the classification of MS potentials for the 
Darboux space of type III given in [9]. At this point, two different proposals for the 
definition of an AD QMS anharmonic oscillator perturbation on this space arise in a 
natural way. The first one consists in the same type of generalization proposed in the 
constant curvature cases described in section 2: 

a + J_ a + J_ ^ \a + J- 

k=l 

2 2 00 / 2 \ 

a + q 2 a + q 2 f-J \a + q 2 J 

A second (and perhaps more natural) possibility arises if we realize that in (3.3) the 
intrinsic oscillator on this curved space is just the "Euclidean" oscillator divided by the 
conformal factor of the metric, which carries all the information concerning the non- 
constant curvature of the space (see [19]). From this perspective, the curved anharmonic 
oscillator would be defined as 

H _ J + + ^ 2 J + g=i ^J- k+1 _ p 2 + ^ 2 q 2 + E£i (37) 

a + J- a + q 2 

In any case, both Hamiltonians (3.6) and (3.7) are QMS and, like the rest of the systems 
presented in this paper, they do have the same set of universal integrals (1.4) coming from 
their sl(2,R) symmetry. 
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